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THEORETICAL  AND  NUMERICAL  SIMULATION  STUDIES 
OF  MIDLATITUDE  F REGION  IRREGULARITIES 

I.  INTRODUCTION 

Temperate  — latitude  spread  F has  been  observed  for  many 
years  (see  Herman,  1966  for  a comprehensive  review!.  Recently, 
a model  which  attempts  to  describe  the  phenomenon  has  been  pro- 
vided by  Perkins  [1973],  Perkins  derived  a set  of  magnetic  field 
line  integrated  moment  equations  (thus  allowing  variations  only 
transverse  to  the  geomagnetic  field;  see  Section  IH  for  the  dynamical 
evolution  of  the  nighttime  F region.  He  showed  that  the  equilibrium 
in  which  the  nighttime  F region  is  supported  by  E x B drifts  is 
unstable,  if  in  addition  to  the  supporting  eastward  electric  field, 
a north-south  electric  field  exists.  The  instability  takes  the  form 

of  rising  and  falling  sheets  of  ionization  and  his  linear  analysis 

-4  -1 

showed  a growth  rate  y « 3 x 10  u.  (where  is  the  ion-neutral 

in  in 

collision  frequency  at  the  F region  peak!.  This  can  provide  growth 
times  of  the  order  of  a few  minutes  to  an  hour  depending  on  whether 
the  F layer  is  high  or  low.  Perkins  showed  that  the  results  of  this 
instability  generation  mechanism  were  able  to  predict,  albeit  based 
on  a linear  analysis,  several  of  the  features  observed  in  spread  F. 

In  Perkins'  derivation  of  the  linear  growth  rate  an  equil- 
ibrium was  taken  where  both  the  field  line  integrated  Pedersen 
conductivity  and  field  line  integrated  ion  number  density  were 
assumed  to  have  no  horizontal  gradients.  However,  as  noted  by 
Perkins  [ 1973],  this  spread  F instability  will  cause  Pedersen  con- 
ductivity gradients  transverse  to  the  magnetic  field  to  be  formed. 
These  gradients  should  then  be  subject  to  the  usual  E x B instability 
(Simon.  1963,  1970;  Linson  and  Workman,  1970!, 

Note:  Manuscript  submitted  February  14,  1975. 
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It  is  the  purpose  of  the  present  paper  to  extend  the  work  of 
Perkins  to  the  nonlinear  development  of  a spatially  varying  ionospheric 
equilibrium.  In  particular,  numerical  simulation  of  Perkins'  [1973] 
moment  equations,  with  an  initial  equilibrium  which  has  a long  wave- 
length variation  of  the  integrated  Pedersen  conductivity  in  the  east- 
west  direction  and  a shorter  wavelength  perturbation  of  this  conductivity 
in  the  north-south  direction,  evolves  to  a quasi-final  nonlinear  state 

with  the  squared  amplitude  of  the  Pedersen  conductivity  fluctuations  vary- 
_2 

ing  as  k (where  k is  the  magnitude  of  the  wavevector).  Such  a power 
law  for  the  F region  irregularities  is  in  good  agreement  with  observa- 
tions (Dyson  et  al. , 1974;  Rufenach,  1972). 

This  paper  is  divided  into  four  sections.  In  Section  II  we 
present  and  discuss  a set  of  moment  equations  first  derived  by  Perkins 
1 1973]  which  describes  the  dynamical  behavior  of  the  nighttime  F region. 
In  addition,  a plausible  equilibrium  state  is  described  and  a linear 
perturbation  analysis  is  carried  out  under  conditions  in  which  F region 
irregularities  can  develop.  A delineation  of  the  differences  between 
our  equilibrium  and  perturbation  analysis  and  that  of  Perkins  [1973] 
is  presented  in  this  section.  In  Section  III  numerical  simulations  of  the 
nonlinear  evolution  of  the  dynamical  equations  for  a portion  of  the 
nighttime  F region  are  presented  and  discussed.  Section  IV  contains 
a summary  and  further  discussion  of  the  results. 


2 


II.  MOMENT  EQUATIONS;  EQUILIBRIUM  AND 
LINEAR  ANALYSIS 


The  derivation  of  the  moment  equations  presented  here  was 
first  carried  out  by  Perkins  [1973].  We  will  not,  therefore,  rederive 
this  set  of  equations,  but  reiterate  the  necessary  assumptions  made 
in  their  derivation.  They  represent  a set  of  moment  equations,  in  the 
electrostatic  limit,  for  the  magnetic  field  line  integrated  ion  number 
density  and  Pedersen  conductivity  of  the  F region.  The  field  line 
integrated  potential  is  derived  self- consistently  from  the  vanishing 
of  the  divergence  of  the  perpendicular  current. 

The  assumptions  made  in  the  derivation  of  the  set  of  equa- 
tions are  as  follows;  (U  the  magnetic  field  is  uniform  and  makes  a 
dip  angle  D with  the  horizontal  direction  (see  Figure  1);  (2)  E region 
contributions  to  the  ionospheric  conductivity  are  ignored;  (31  recom- 
bination and  ionization  are  neglected;  (4)  the  neutral  atmosphere  con- 
sists of  a single  species  distributed  vertically  with  a scale  height  H ; 
(5)  both  the  falling  of  ions  and  horizontal  drifts  due  to  gravity  are 
included;  (6)  the  plasma  is  high  enough  so  that  the  ion  cyclotron  fre- 
quency is  greater  than  the  ion-neutral  collision  frequency  • 

in’ 

(7)  the  ions  and  electrons  are  isothermal;  and  (8)  electric  field  and 
neutral  winds  give  rise  to  Pedersen  currents  which  can  support  the 
ionosphere  by  j^x  B forces.  The  geometry  for  the  problem  is  given 
in  Figure  1. 

The  equations  can  now  be  written  as; 


n 


■E 


c cos  D 
B H 

n 


(1) 

(2) 


3 


(3) 


V • (E  V Ecos  D (■ 

I I ^ 


2T 


J.  J. 


eH 

n 


ge  cos  D 

n 
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where  equations  (1)  and  (2)  represent  the  transport  equations  for  the 
integrated  number  density  and  Pedersen  conductivity  respectively. 
The  symbols  used  in  equations  (li.  (2').  and  (31  are  defined  as  follows; 

N = field  line  integrated  ion  number  density 
E = field  line  integrated  ion  Pedersen  conductivity 
= field  line  integrated  potential 
g = gravitational  acceleration 
B = ambient  geomagnetic  field 
n = ion  cyclotron  frequency 
= neutral  scale  height 
m = ion  mass 
D = magnetic  dip  angle 
e = electron  charge 
T = temperature  in  units  of  energy 
J = direction  transverse  to  B. 

Equation  (3)  is  the  equation  for  the  potential  obtained  from 


where  J is  the  field  line  integrated  current.  It  should  be  emphasized 

-L 

that  equations  (1)  - (3)  are  two-dimensional,  i.e.,  represent  variations 
transverse  to  the  ambient  geomagnetic  field. 


west  direction  will  be  formed  and  it  must  be  shown  that  these  gradients 
persist  throughout  the  night.  To  show  that  these  gradients  do  exist  we 
will  show  that  the  set  of  moment  equations  yield  a self-consistent 


V • J = 0 
-L  J. 


(41 


In  the  predawn  and  post-sunset  hours,  gradients  in  the  east- 
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equilibrium  when  gradients  of  N and  exist  in  the  east -west  or  y 
direction.  This  represents  a physical  situation  which  is  favorable 
for  the  development  of  the  ^ x B gradient  drift  instability  provided 
the  east-west  gradients  persist  and  the  ion-neutral  relative  drift 
is  greater  than  a particular  threshold  value. 


To  do  this  assume  that 


— — = - — = 0 and  — — = -E  , 
dx  dx  dx  ox 


where  subscript  zero  denotes  equilibrium  quantities.  Thus,  there 
are  no  north-south  gradients  in  N and  Z).  and  a constant  north-south 
electric  field  exists.  Under  these  conditions  the  set  of  moment 
equations  yield 


, ^^o  / gcosD  ^^ox\_  ecgsin^D  , v*  ccosD 

di  dy\n~B  llBH  ^o  ^ -^o  dy  BH 

n n 


5N  5N  / 

0^  ^ o /g  cosD 

^ dy  \ n 


c E \ 
ox  \ 

B / 
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(6) 


d , ^‘^o  ) 2 T *^"^0  ^^o  gecosD 

^y  '“^o  ■JjT  e dy^  ~ 


(7) 


The  equilibrium  that  we  seek  exists  in  the  frame  drifting  with  the 
ions  in  the  y direction,  therefore,  we  transform  to  a frame  drifting 
with  a velocity 

c E 

_ g cos  D ox  2T  c cos  D 

y n ■ B ^ e B H 


Equations  (5).  (61.  and  (7)  become 

al3o  0^0  I 2T  ccosD\_  ecgsin^D  ,,  v*  ccosD 

dt  ^ dy  \ e BH  I UBH  ^o  ^o  &y  BH 
' n ' n n 
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(10) 


/_  2;^  ccosD\ 

5 1 5y  y e B j 
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d ^ V ^ 2T  ^ g e cos  D 
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(11) 


The  first  integral  of  equation  (11)  yields 


9^0 

9y 
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N 


o g e cos  D 


ft 


(12) 


which  is  actually  a condition  on  the  east-west  component  of  the  current, 
involving  a Pedersen  conductivity  gradient  in  the  east-west  direction. 

At  this  point  we  assume  that  the  gradient  of  N is  very  small. 
Thus,  the  gradient  of  Pedersen  conductivity  is  due  largely  to  a shift 
in  the  altitude  along  a field  line  of  the  F region  ion  density  peak  as 
one  proceeds  in  an  east-west  direction.  The  integrated  Pedersen 
conductivity  increases  as  the  ions  on  a field  line  move  down  the  field 
line  to  regions  where  the  ion-neutral  collision  frequency  increases 
because  the  number  of  neutrals  increase  with  lower  altitude.  Using 
the  assumption  of  small  gradients  in  N and  substituting  equation  (12) 
into  equation  (9)  yields 


aDo 

at 

aN 

o 

at 

Setting 


2 2 

e c g sin  D , c e g cos  D „ ^ c cos  D 

ftBH  ^o  ftBH  ^o"  BH 
n n n 


0 


at 


0 yields  a value  for  the  integration  constant  K. 


K 


e g N 
^ o 

ft  cos  D 


(13) 

(14) 


(15) 
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Now  that  it  has  been  shown  that  an  equilibrium  exists  when 
there  is  a gradient  of  the  Pedersen  conductivity  in  the  east -west 
direction,  it  remains  to  consider  localized  perturbations  about  that 
equilibrium.  These  perturbations  are  of  the  form  [e^^^x  ^ 

where  k^,  and  the  linearized  equations  are  given  as 

^ ^'(.2T  ccosD,  aEo  ( ^ c . 

di  dy  ^ e B ^ ay  ax  ' '^y%B^ 


c cos  D (90'  c cos  D 

- Z-  a B H 


9y 


n 


“oay  B H 


n 


*’)  + ^ ^ + Jf-)  vfi;' 

n 


(16) 


= 0 


(17) 


where  the  primes  denote  perturbed  quantities. 

Substituting  the  functional  form  of  the  perturbations  into 
equations  (16)  and  (17)  yields 


- i k 


_ 2T 
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5 00 
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• 1 f 2T 
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(19) 


7 


Substituting  equation  (19)  into  (18)  and  assuming  w = + iy  yields 

the  expression  for  the  growth  rate 
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c cos  D 
ay  B H 
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‘'x  c 
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c cos  D 
y B H 
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\ dy 


4-  k 
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5(f>0 

5x 


k cos 

X 
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ay  e ' X y 


(20) 


The  above  expression  (20)  represents  the  growth  rate  for  a 
mode  which  in  general  lies  at  an  arbitrary  angle  to  the  east  west 
direction.  To  gain  a better  physical  insight  into  the  phenomena  we 
will  restrict  our  attention  to  those  modes  whose  wave  vectors  lie 
along  the  x or  north-south  direction  (see  Figure  l.b).  Thus  “ 0 
and 


y 


£ / Bg  cosD  _ \ 1 5i:o 

B \ lie  ■ ox/^^  dy 


COS  D 


z 

Eo 


(21) 


where 

.2^  N 
g ^ e g sin  D o 

^oy  ft  cos  D 

and  it  is  assumed  that  the  functional  form  of  is 

o 


E 


(22) 


Eq  = E + L^(y),  ^ 
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where  Y,  is  the  spatially  averaged  part  of  the  equilibrium  Pedersen 
conductivity  and  contains  the  east-west  spatial  variation,  the 
condition  that  a perturbation  grows  becomes 


c / gcosDB 
B\  n c 


1 *^^o 


cos  D E 


(23) 


Thus,  it  has  been  shown  that  an  equilibrium  can  exist  when  an 
east-west  gradient  in  the  Pedersen  conductivity  is  present  and  that 
perturbations  about  that  equilibrium  with  wavevectors  in  the  north- 
south  direction  will  grow.  This  equilibrium  can  be  compared  with 
the  results  of  Perkins  [1973].  Perkins  found  that  if  one  chooses  an 
equilibrium  of  the  form 


= const 

N = const 
o 

E = const 

■ 2 ^ N 
_ e g sin  D o 

^oy  " ncos  D 


and  perturbs  about  that  equilibrium,  those  modes  with  wavevectors 
lying  between  the  zeroth  order  electric  field  and  the  east-west  direction 
would  be  unstable.  The  marginally  stable  mode  has  a wavevector 
parallel  to  the  east -west  direction. 
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III.  NUMERICAL  SIMULATION  RESULTS 


The  simulations  were  initialized  with  perturbations  in  the 
north-south  direction,  as  follows: 


L(x,y)  = + ^’(x); 


(24) 


= E[1  + sin  k^y], 


y = 


2tt 


30  km 


(25) 


^'(x)  = Ee  ^ sin  k x, 


X = 


2tt 


10  km 


(26) 


oy 


2 

eg  sin  D 
11  cos  D 


N 


(27) 


E = tan  0 E , 
ox  oy 


(28) 


with  9 given  below.  The  moment  equations  can  be  written,  after  trans- 
forming to  a frame  moving  with  velocity 


oy  ^ 
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(E  c 

( B-- 


gcosD  V A 

fi 


(29) 
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and  defining 


Y.=  X)(x,  y)  - Z)  _ ^ 

L ^ 


(30) 


V<f}  = -E  + V<^' 


(31) 


as: 
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5E 


V^'  <f>'  + ^ V </>*  + ^ ^ cos 


7^  ,ht 


1 .E 


i+E 


d( 

\ e H e / 


(32) 


^ _E  .5^  , 0 


i+E 
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(33) 


The  system  of  equations  (32)  and  (33)  for  E and  (/>  was  solved 
numerically  on  an  IBM  360/195  computer.  The  transport  equation  was 
solved  using  a two-dimensional  Lax-Wendroff  scheme,  while  the  poten- 
tial was  obtained  using  a modified  ADI  algorithm.  A two-dimensional 
rectangular  32  x 32  grid  with  mesh  size  Ax  = Ay  = 1 km  was  chosen. 

The  X and  y axes  corresponded  to  the  north-south  and  east-west  direc- 
tion, respectively.  The  boundary  conditions  were  periodic  in  both 
directions.  The  integrated  ion  number  density  was  constant  throughout 


11 


the  various  simulations  (c,f.  Perkins,  1973).  The  parameters 
chosen  for  the  simulation  were:  D = 45°; 

N = 5 X 10^^/. 45  = 1/9  X 10^^  (mhos  - cm^)"\  resulting  in 

® ® -4 

^in  /fl  = 2.8  X 10  ; B = 0.5  gauss;  m = 16  m (m  is  the  proton 

o P P 

mass);  T = 916  K;  = 55  km.  These  parameters  were  chosen  to 

correspond  to  an  altitude  -350  km  (see  Johnson,  1961,  and  Rishbeth 

and  Garriott,  1969). 

The  results  of  this  investigation  are  best  described  within 
the  context  of  Figures  2 through  5.  Figures  2 and  3 show  the  evolu- 
tion in  time  of  an  F region  slab  for  two  cases  in  which  the  ratio  of 
north-south  to  east-west  electric  field  was  changed.  Figure  2 repre- 
sents the  damped  solution  f or  the  case  for  which 
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= cos  D 


max 
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where  E 


oy 


eg  sin  D 
Rcos  D 


^ and  we  took  E =4.1  mV/m  and 

o - . ^ 

9 = 61  (Eqj^  =1.7  mV/m).  The  linear  growth  rate  for  this  case 


(where 
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aZo  is  a maximum)  is  given  as: 
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This  situation  yields  growth  only  in  regions  where 


contours  seen  at  the  extreme  right  hand  side  of  the  1450  second  plot 


of  Figure  2. 

Figure  3 represents  a physical  situation  in  which  the  north- 
south  mode  is  a growing  mode.  For  this  case 


Since  » g cos  D m/e  the  linear  growth  rate  for  this  case  is  given 
approximately  as: 


Therefore,  the  mode  only  grows  in  regions  where  the  gradient  in 
Pedersen  conductivity  is  negative.  This  is  born  out  in  Figure  3.b 
where  growth  is  exhibited  in  the  central  region,  i.e.,  where  the 
initial  sinusoidal  y variation  is  favorable. 


scale  structures  at  late  times  from  the  initial  large  scale  structures 
(Figure  3c).  This  shows  that  there  is  a cascading  of  energy  from 
long  wavelength  modes  to  shorter  wavelength  modes.  This  cascading 
of  energy  and  an  attendant  power  law  for  the  power  spectrum  of  the 
fluctuations  is  reminiscent  of  the  recent  computer  simulation  results 


max 


Here  E =4.1  mv/ m and  E =21.6  mV/ m {9  =11°). 
oy  ox 


y = 


E 


o 


An  important  point  to  note  is  the  development  of  the  smaller 
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for  type  II  irregularities  in  the  equatorial  E region  electrojet  given  by 

McDonald  et  al.  [1974].  In  both  the  E region  studies  and  the  present  F 

region  studies  the  basic  generating  mechanism  is  a form  of  the  E x B 

gradient  drift  instability  (Simon.  1963;  1970;  Linson  and  Workman, 

1970).  However,  in  the  two  cases  the  dynamical  equations  governing 

each  region  are  different.  Figure  4 shows  the  development  of 

(823/S)j^  for  three  different  modes  in  the  system.  The  modes  labeled 

0 and  1 are  the  east-west  and  north-south  modes,  respectively,  with 

which  the  system  was  initialized.  We  note  that  mode  1 saturates  at 

-4.8%  of  the  ambient  level.  The  mode  labeled  2,  however,  is  the 

first  harmonic  of  mode  1.  Mode  2 is  initially  lost  in  the  background 

noise,  but  eventually  appears  and  grows,  indicating  that  smaller  scale 

structures  are  developing.  Many  other  shorter  wavelength  modes 

develop  as  the  simulation  proceeds;  however,  we  chose  only  to  display 

the  first  harmonic  of  mode  1.  Figure  5 represents  a direction  averaged 

power  spectrum  for  the  Pedersen  conductivity  contours  depicted  in 

Figure  3c  (t  = 669  seconds).  A least-squares  fit  of  the  numerical 

_2 

simulation  results  (dots)  with  a power  law  curve  yielded  a k spectrum. 
The  falling  off  of  the  curve  of  the  short  wavelengths  is  due  to  the  num- 
ber of  grid  points  used  (32x32  km  with  Ax  = Ay  = 1 km),  i.e.,  to  com- 
puter resolution  for  this  particular  simulation.  This  result  is  in  good 
agreement  with  the  power  spectrum  obtained  by  in  situ  measurements 
from  OGO  6 of  ion  number  density  fluctuations  made  by  Dyson  et  al. 
[1974].  These  authors  reported  that  the  square  of  the  amplitude 
fluctuations  obeyed  a power  law  of  the  form  k"*^  with  n close  to  2 and 
an  average  value  of  1.9  for  wavelengths  between  70m  and  7 km. 

The  power  law  result  from  the  computer  simulation  is  also  in 
agreement  with  the  power  spectrum  results  deduced  by  Rufenach  [1972] 
from  ionospheric  scintillation  observations  (at  Boulder,  Colorado) 
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from  the  radio  source  Cygnus  A analyzed  at  26  MHz.  Rufenach  showed 
that  the  ix)wer  spectrum  for  the  irregularities  in  the  wavelength  regime 
from  0.6  to  4 km  was  ock  where  n typically  varied  between  2 and  3. 
The  total  for  t = 669  seconds  is  ~10%  which  is  not  at  odds 

with  the  in  situ  measurements  of  Dyson  et  al.  [1974]. 

A phenomenon  associated  with  spread  F is  ducting  of  radio 
waves  along  field  lines.  The  presence  of  localized  enhancements  and 
depletions  in  the  background  Pedersen  conductivity  seen  in  Figure  3 
represents  a physical  situation  favorable  to  ducting  of  e.m.  waves. 
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IV.  SUMMARY  AND  DISCUSSION 


We  have  used  the  nighttime  temperate  latitude  F region  equa- 
tions for  field  line  integrated  Pedersen  conductivity  and  potential 
developed  by  Perkins  [ 1973]  to  extend  his  work  to  the  nonlinear  develop- 
ment of  a spatially  varying  ionospheric  equilibrium.  In  particular,  we 
have  used  a long  wavelength  variation  of  Pedersen  conductivity  in  the 
east-west  direction  as  part  of  the  equilibrium  to  which  we  add  at 
t = 0 a shorter  wavelength  north-south  perturbation  of  conductivity. 

Such  a system,  under  appropriate  conditions  [see  equation  (23)],  is 
subject  to  the  usual  E x B gradient  drift  instability.  The  nonlinear 
evolution  of  such  a system  produces  a power  spectrum  for  the  Pedersen 
conductivity  fluctuations  which  obeys  the  relation 
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This  nonlinear  distribution  of  energy  involves  a cascade  process  in 
which  energy  is  transferred  from  longer  wavelength  modes  to  shorter 
wavelength  modes.  This  result  is  in  agreement  with  in  situ  measure- 
ments of  F region  irregularities  Dyson  et  al.  [1974],  with  observations 
from  scintillation  measurements  taken  at  Boulder,  Colorado  [Rufenach, 
1972],  and  with  incoherent  scatter  measurements  from  Arecibo 
[Rino  and  Livingston,  1974]. 

The  initial  orientation  of  our  long  wavelength  equilibrium  and 
perturbation  produces  structures,  in  the  nonlinear  state,  which  are 
sheet-like  and  more  aligned  in  the  east-west  direction,  i.e.,  are 
axially  asymmetric  in  the  plane  perpendicular  to  the  ambient  geomagnetic 
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field.  Such  east- west  structures  have  been  observed  by  Singleton 
[1970]  during  the  winter,  at  midlatitudes,  by  using  scintillation  data 
from  the  Explorer  22  beacon  experiment.  However,  at  this  point 
the  shape  of  the  F region  irregularities,  determined  by  scintillation 
observations,  is  not  entirely  clear.  A good  portion  of  the  time  the 
irregularities  appear  to  be  axially  symmetric  (Singleton,  1970; 
Moorcroft  and  Arima,  1972).  Moorcroft  and  Arima  [1972]  also  give 
evidence  for  elongation  more  in  the  north-south  direction  as  does 
Singleton  [ 1970]  for  seasons  other  than  winter  (it  should  be  pointed 
out  that  Moorcroft  and  Arima’s  data  was  not  taken  in  winter  and  so 
is  not  inconsistent  with  Singleton’s  findings).  Our  computer  simula- 
tion results.  Figure  3,  would  depict  an  average  k for  the  final  state 
more  in  the  north-south  direction.  Phase  space  plots  in  k space  sup- 
port this  idea,  i.e.,  that  the  final  power  spectrum  has  a preferential 
direction  (north-south).  This  is  unlike  the  numerical  simulation 
results  of  the  type  II  equatorial  E region  irregularities  which  show 
that  the  k vectors  are  isotropized  in  the  final  state  (see  McDonald 
et  al.  1974b).  This  is  probably  due  to  Hall  current  effects  which 
operate  effectively  in  the  E region  but  not  in  the  F region.  For 
the  initial  conditions  chosen  for  our  simulation,  this  implies  more 
of  a one-dimensional  steepening  process  rather  than  a true  two- 
dimensional  turbulence  process  (as  is  the  case  of  type  II  equatorial 
E region  irregularities;  see  McDonald  et  al.  1974a)  for  the  result- 
ing nonlinear  quasi -final  state.  We  believe  that  if  we  choose  a 
long  wavelength  north-south  equilibrium  and  a perturbation  in  the 
east-west  direction,  the  resulting  nonlinear  development  of  the 
E X B instability  would  show  a north-south  preference.  We  are 
continuing  further  studies  along  these  lines.  Moorcroft  and  Arima 
[1972]  have  pointed  out,  in  a qualitative  fashion,  that  the  E x B 
gradient  drift  instability  is  a good  candidate  for  the  axial 
asymmetry. 
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We  do  not  wish  to  imply  that  our  results  can  explain  all 
spread  F observations.  Rather,  we  have  shown  how  the  ^ x B 
instability,  using  Perkins'[  1973]  F region  equations,  can  produce  F 
region  irregularities  which  are  in  agreement  with  certain  aspects  of 
the  observations  of  these  irregularities.  Perkins'  instability  mechanism 
had  a growth  rate  proportional  to  (not  an  B instability,  since 

his  initial  equilibrium  had  a constant  value  for  the  F region  Pedersen 
conductivity).  He  compared  his  results  with  experimental  observations 
of  spread  F (see  Herman,  1966  for  a comprehensive  review  of  spread 
F).  He  concluded,  based  on  linear  theory,  that  the  mechanism  was  in 
agreement  with  spread  F observations.  Perkins'  instability  mechan- 
ism required  a north-south  electric  field,  as  does  our  mechanism.  His 
instability,  being  inversely  proportional  to  the  F region  ion-neutral 
collision  frequency,  favored  instability  when  the  F layer  was  high. 

Our  instability  criterion,  equation  (23),  while  involving  the  north-south 
and  the  east-west  electric  field,  also  involves  the  neutral  scale  height, 
Hj^  (this  comes  from  Perkins'  model  equations  where  the  ion-neutral 
collision  frequency  is  taken  as  having  exponential  dependence  with 
this  scale  height).  The  condition  shows  that  for  a fixed  ratio  of  elec- 
tric fields  the  instability  is  stronger  for  larger  values  of  H^,  i.e., 
increasing  F region  altitudes.  Also  of  interest  is  the  fact  that 
neither  instability  mechanism  makes  use  of  a vertical  gradient  in  ion 
density. 

In  concluding,  we  have  shown  by  a numerical  simulation  fol- 
lowing Perkins'  [1973]  nighttime  F region  equations  in  time  into  the 
nonlinear  regime,  that  F region  irregularities  can  be  produced  via  the 
usual  E X B instability.  However,  we  believe  that  further  work  on  the 
generation  and  nonlinear  development  of  F region  irregularities  is 
warranted  in  order  to  gain  an  improved  understanding  of  the  relation- 
ship between  these  irregularities,  spread  F and  scintillation  phenomena. 
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Fig.  4 — Pedersen  conductivity  fluctuation  amplitude  squared  for  several 
modes  as  a function  of  time  for  the  case  depicted  in  Fig.  3.  The  main 
initial  east-west  long  wavelength  spatial  variation  is  labeled  0 (squares)  and 
the  initial  shorter  wavelength  north-south  perturbation  is  labeled  1 
(dots).  The  label  2 (triangles)  is  the  second  harmonic  of  1. 
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Fig.  5 — Direction-averaged  power  spectrum  of  the  Pedersen  conductivity 
fluctuations  at  669  seconds  depicted  in  Fig.  3,  plotted  on  a log-log  scale. 
On  the  wavenumber  scale,  ko  is  the  initial  long  wavelength  east-west  vari- 
ation. The  power  law  k~2  was  obtained  by  a weighted  least  squares  fit 
to  the  simulation  results  (dots). 
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